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In Affleck-Dine baryogenesis, the observed baryon asymmetry of the Universe is generated through 
the evolution of the vacuum expectation value (VEV) of a scalar condensate. This scalar condensate 
generically fragments into non-topological solitons (Q-balls). If they are sufficiently long-lived, they 
lead to an early matter domination epoch, which enhances the primordial gravitational wave signal 
for modes that enter the horizon during this epoch. The sudden decay of the Q-balls results in a 
rapid transition from matter to radiation domination, producing a sharp peak in the gravitational 
wave power spectrum. Avoiding the gravitino over-abundance problem favours scenarios where the 
peak frequency of the resonance is within the range of the Einstein Telescope and/or DECIGO. This 
observable signal provides a mechanism to test Affleck-Dine baryogenesis. 


The asymmetry between matter and anti-matter is a 
cornerstone puzzle of modern particle cosmology, as the 
Standard Model fails to provide an explanation [IH3]. 
An elegant paradigm for explaining the asymmetry is 
the Affleck-Dine mechanism [3}6]. Supersymmetric the- 
ories generically have flat directions [5] [7], which have 
non-zero baryon or lepton number. During inflation, 
a scalar condensate generically develops in these direc- 
tions, whose non-zero vacuum expectation value (VEV) 
spontaneously breaks C and CP. At the end of infla- 
tion, a baryon and/or lepton asymmetry is generated 
as the VEV coherently evolves and the condensate frag- 
ments [8]. These resulting clumps may be long-lived non- 
topological solitons (Q-balls) [9112], carrying either lep- 
ton or baryon number [13]. This global charge is trans- 
ferred to Standard Model particles when the Q-balls de- 
cay. 

However, the Affleck-Dine mechanism is generically 
a high-scale phenomenon which does not require a low 
SUSY-breaking scale, making it difficult to confirm ob- 
servationally. In this letter, we argue that a broad class 
of Affleck-Dine models significantly enhance the primor- 
dial gravitational wave power spectrum. This provides a 
novel mechanism to test or constrain these models. 

Generically, Q-balls produced through the fragmenta- 
tion of the Affleck-Dine condensate are large and long- 
lived. Consequently, they evolve as non-relativistic mat- 
ter, and may eventually come to dominate the energy 
density of the Universe. If the Q-balls decay rapidly, 
there is a sudden change in the equation of state for 
the Universe. This results in rapidly oscillating scalar 
perturbations, which enhances the primordial gravita- 
tional wave spectrum from inflation. This is analogous 
to the poltergeist mechanism, in which the sudden de- 


graham.white@ipmu.jp 
t Ipearce@psu.edu 


usenko@g.ucla.edu 


Author order determined by coinflip 


q 


cay of black holes also enhances the gravitational wave 
spectrum [14]. This is in contrast to the production of 
gravitational waves during the fragmentation of the con- 
densate, as typically the condensate is a small fraction 
of the initial total energy [16]. Our proposed test is 
also complimentary to constraints on the ratio of scalar 
to isocurvature perturbations from D-term inflation [I7] 
and those from the backreaction of the Affleck-Dine po- 
tential on the inflaton potential, which can impact the 
cosmic microwave background [18]. 

In this letter, we first argue that Affleck-Dine scenarios 
generically have this epoch of early matter domination, 
and secondly, the Q-ball decay rate is sufficiently fast to 
enhance the gravitational wave spectrum. In particular, 
the sudden transition avoids the suppression that occurs 
in a gradual transition like Moduli decay [19]. Analytical 
arguments and simulations show that the Q-ball mass 
distribution is sharply peaked [8] 20H22]. Secondly, the 
charge quanta inside the Q-balls decay to fermions. Q- 
balls decay when the decay rate per unit charge is larger 
than the Hubble parameter. We show below that this 
is suppressed by a surface area to volume factor, and 
therefore the decay rate per unit charge accelerates as 
the Q-ball decays, similar to black hole decay. 

Furthermore, avoiding an overabundance of gravitinos 
results in a gravitational wave spectrum that is at suffi- 
ciently low frequencies to be observed. Finally, although 
in this work we make no statistical claims, we present 
several points in parameter space where the Q-balls are 
sufficiently long lived to dominate the energy density and 
produce a detectable gravitational wave signal. Thus, 
if such a signal is observed, we can narrow the cause 
down to two known scenarios - an early period of Q-ball 
domination, which is likely a consequence of Affleck-Dine 
baryogenesis, or an early period of light primordial black 
hole domination [14] [23] [24]. 

Q-ball Induced Early Matter Domination: dur- 
ing inflation, the field ® acquires a vacuum expectation 
value when averaged over super-horizon scales [4 [25}{28). 
At the end of inflation, it relaxes towards its equilibrium 
value as the field fragments to form Q-balls [3] [8]. Dur- 


ing the relaxation process, a charge excess is produced 
as the field VEV follows a curving path in field space, 
which is biased either clockwise or counter-clockwise by 
the small CP-violating operator. However, as a higher 
dimensional operator, it will also be sensitive to the ini- 
tial post-inflationary VEV, which is subject to random 
fluctuations during inflation. Consequently, some Hub- 
ble patches will have an excess of Q charge while other 
have an excess of Q charge. Therefore, there are sym- 
metric and asymmetric components to the initial Q-ball 
density. 

After fragmentation, most of the condensate’s initial 
energy is contained in Q-balls rather than individual par- 
ticles, particularly if the couplings between the scalar 
field and fermions is small [8]|20]. If the asymmetric com- 
ponent is small (as is expected due to the smallness of the 
observed baryon asymmetry), the symmetric component 
must then be large. We parameterize the asymmetric 
component of the Q-ball energy density 


non 
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and we expect r to be within an order of magnitude of 
the baryon asymmetry. (This can also be understood as 
a consequence of a highly elliptical orbit during relax- 
ation, which simulations connect to a large symmetric 
component [22].) 

In the thin wall regime, the vacuum expectation value 
inside the Q-ball can be found by minimizing V (7) /®? 
where ® parameterizes the flat direction in the Affleck- 
Dine potential. The energy per unit charge of a single 
Q-ball is given by 


> (2) 


where v is the VEV inside the Q-ball (we discuss specific 
potentials in Appendix B) The total initial energy in 
Q-balls after fragmentation is 


PQ = Qowno (3) 


where no ~ NQH§ is the initial number density of the 
Q-balls and Qo is their initial charge. Simulations sug- 
gest Ng ~ 1000 for gravity-mediated SUSY scenarios 
and Ng ~ O(1) for most gauge-mediated scenarios, if 
higher dimensional operators are negligible. However, in 
this scenario the resulting Q-balls are in the thick wall 
(as opposed to thin wall) regime [27]. Although we focus 
on thin wall Q-balls in this work, we note that scaling ar- 
guments suggest thick wall Q-balls are longer lived and 
thus also can induce early matter dominated epochs (see 
the supplementary material [1}. 

From this, it is straightforward to derive an expression 
for the initial charge of the Q-balls in terms of the initial 
baryon asymmetry Ypo, r, and the reheating temperature 
To, 
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Because the total charge in each Hubble volume is dis- 
tributed into a small number of Q-balls during fragmen- 
tation, the initial Q-balls tend to have large charge; in 
our benchmark scenarios, the initial charges are above 
107°. Consequently, they will travel at non-relativistic 
speeds in the post-inflationary plasma. Then, if the Q- 
balls live long enough, they dominate the energy density 
of the universe. We can approximate the temperature of 
matter-radiation equality in the limit where Q-ball decay 
is negligible as 
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Teg ~ (5) 


Although long-lived, the Q-balls produced by the frag- 
mentation of the Affleck-Dine condensate are not sta- 
ble since their conserved charge must be transferred to 
Standard Model particles. The sfermions carrying the 
charge can decay to a quark (or lepton) and neutralino 
or chargino. Expressions for the relevant coupling can 
be found in Ref. [29], although we will parameterize the 
vertex in terms of an effective Yukawa coupling Yer. 


This decay happens only at the surface of the Q-ball, 
for one of two reasons. First, if the VEV of the squark or 
slepton fields inside the Q-ball (which carry the charge) 
is significantly larger then the energy per unit charge w, 
then the large induced fermion masses forbid the decay 
inside the Q-ball. The induced masses of the Standard 
Model fermions have magnitude gv, where g = g3 for 
quarks if the Q-ball is made of squarks and g = ge for 
leptons if the Q-balls is made of sleptons [30]. Therefore, 
if gu = w, the decay occurs only at the surface of the 
Q-ball, where the VEV drops to zero. 


Alternatively, if the decay is not forbidden, decays in 
the interior of the Q-ball rapidly fill up the Fermi sea. 
Thereafter, the Q-ball quanta decay only at the surface 
as long as the diffusion time, tp ~ 3R?/A, is sufficiently 
long. The mean free path is À ~ 1/oygn, where number 
density n = 3Qo/(47R?) refers to the density of scalar 
quanta inside the Q-ball. The diffusion time is shortest 
for the highest momentum, which is ~ w/2 when the 
decay is energetically forbidden. The diffusion time can 
then be approximated using the scattering cross section 
Os ~ Gi /(w/2)?, where gi € (gy, 92,93) depending on 
the Standard Model fermion and sfermion involved. For 
the benchmark points presented below, decays inside the 
Q-ball are suppressed for the first reason. 


Regardless of the reason, the Q-ball evaporation rate 
is suppressed by the ratio of the surface area to the vol- 
ume. Because the radius scales as Q!/°, the decay rate 
then scales as Q?/* instead of as Q. A Q-ball decays once 
the decay rate per unit charge is larger than the Hub- 
ble parameter. When decays occur only on the surface, 
Te-pan/Q x Q713, which means that it accelerates as 
the Q-ball shrinks. Therefore, Q-ball decay is an effec- 
tively instantaneous decay, similar to black holes [74] [BI]. 


The charge depletion per unit time per unit area of a 


Q ball obeys the equation [32] 
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where v is the field value of the condensate and yeg is the 
effective Yukawa coupling, accounting for mixing angles. 
In the thin wall limit, 
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which gives 
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For the Q-balls to decay after dominating the energy den- 
sity, we must require [Q_pan/H Q| r-r > Approxi- 
=Toa 


mating the left side at Q = Qo, we find the condition 
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The large symmetric component r ~ Ypo is vital due 
to the Yg 3 factor, which would otherwise make this 
condition difficult to satisfy. As expected, this prefers 
small Yukawa couplings, which result in long-lived Q- 
balls. 

We emphasize that for our numerical analysis, we 
solved the differential equation dQ/dt = —TQ-Ba1n(Q, T). 
We also note that the decay of the Q - balls dilutes the 
initial baryon asymmetry, resulting in a final asymmetry 
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AY, —3/4 
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where Tuec is the temperature at which the Q-balls decay 
and Yp = 8.59 x 1071! as given by Planck [33]. 

Because the Q-ball mass fraction is sharply peaked 
at a single value and the decay is effectively instanta- 
neous, the scale factor and Hubble approximately obey 
step function solutions 
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where ņ is the conformal time; npr is specifically the con- 
formal time at which radiation domination recommences. 

Gravitational waves: we assume inflation generates 
a primordial scalar power spectrum of the form 
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Pe(k) = O (kint — k) As ~ (12) 


for some cutoff scale king, Ns being the spectral tilt, k, 
being the pivot scale and A, being the amplitude at the 


pivot scale. We take typical values of A, = 2.1 x 1079, 
ns = 0.97, ką = 0.05 Mpe~* [83]. 

Scalar perturbations grow with the scale factor dur- 
ing any matter domination epoch, including the Q-ball 
dominated epoch, which can in turn induce gravitational 
waves [31]. Our analysis of the induced gravitational 
wave signal follows [19] and therefore we similarly work 
within the conformal Newtonian gauge and assume Gaus- 
sian curvature perturbations. 

If matter domination is sufficiently long, then pertur- 
bations at small scales can enter the non-linear regime 
where a linear analysis is insufficient. Such non-linearities 
become important at scales kyu ~ 470/nr, where nr is 
the conformal time at which the Q-ball-caused matter 
domination era abruptly ended. In this work, we neglect 
the non-linear regime and therefore restrict ourselves to 
points in parameter space at which the maximum co- 
moving mode enhanced by early matter domination sat- 
isfies kmax < 470/nr. We note that there may still be 
detectable gravitational wave signals in the parameter 
space where this is not satisfied, although we leave the 
analysis of the non-linear regime to future work. 

Using the step function approximations given above, 
the power spectrum of gravitational waves at conformal 
time 7) is 
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where the time averaged power spectrum of the induced 
gravitational waves is related to the scalar (curvature) 
power spectrum by 
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In the above, the time dependence of the gravitational 
waves is 
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a(n) 
(15) 
where the Greens function is the solution of the equation 
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and the source function has the form 
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In these equations w is the equation of state parameter 
and ® is the transfer function of the gravitational poten- 
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FIG. 1. the gravitational wave signal for three benchmark 
scenarios, which have effective Yukawa couplings similar to 
the Standard Model bottom quark (red, dotted), up quark 
(olive, dot-dashed), and electron (black, dashed). These 
clearly produce signals within the reach of future experiments, 
which were taken from Ref. for DECIGO with 3 units 
and an observation time of 1 year, Ref. for LISA with 
an observation time of 4 years, Ref. [40] [41] for THEIA with 
an observation time of 20 years, Ref. for Einstein 
Telescope with an observation time of one year, Ref. [44] for 
the Cosmic Explorer and Ref. for SKA. 


tial, which obeys the evolution equation [35] 
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For a sufficiently quick transition from matter to radia- 
tion domination, we can use the analytic formulae for the 
gravitational wave power spectrum in Ref. [31] which we 
give in the supplementary material [3 

This rapid transition is necessary to produce the sharp 
peak through the “poltergeist” mechanism [14|. Dur- 
ing the early matter domination epoch, density pertur- 
bations in non-relativistic Q-ball modes grow and form 
overdensities. The Q-ball decay, which is rapid as com- 
pared to the Hubble time, converts these overdensities 
into relativistically moving sound waves, which serve as 
sources of gravitational waves. Gravitational waves ex- 
hibit a rapidly growing resonance mode which is ampli- 
fied by interactions with a sound wave comoving at a 
certain angle [14 [BI] [36]. This resonance results in a 
dramatic enhancement at a certain frequency, as can be 
seen in our Fig. It is important that the transition 
to radiation domination is rapid, because otherwise the 
overdensities dissolve gradually and do not result in any 
relativistically moving modes. 

Results: we present the gravitational wave signal for 
three sample points in parameter space in Fig.|1} These 
were chosen to have Yukawa couplings similar in size to 
those in the Standard Model; the precise values of the 
parameters are given in Table [I] To retain generality, we 
specify the VEV v and energy density per charge w of 
the Q-balls, instead of specializing to a particular poten- 


tial. Gravity and gauge-mediation models which produce 
Q-balls with these properties are discussed in the supple- 
mentary material 


Calculated quantities, such as the equality and decay 
temperatures, for these benchmark points are given in 
Table [T] We note that since w is within one order of mag- 
nitude of To, the temperature at which Q-balls are pro- 
duced, it is self-consistent to neglect finite temperature 
corrections to w, which are induced via loop corrections. 
The fragmentation of the Q-ball condensate determines 
the initial charge as given by Eq. (4), which produces the 
scaling Q ~ 10°4(Ygo/r) - (10° GeV /To)?. The resulting 
values are within the typical range (see, e.g., ys 


The observable range is controlled by the proposed fre- 
quency sensitivity of future gravitational wave detectors, 
with higher frequencies probing higher temperatures. At 
present, the highest frequency gravitational wave detec- 
tors with enough sensitivity are the Cosmic Explorer 
and the Einstein Telescope although higher fre- 
quency proposals are a promising work in progress [48]. 
We see that LISA has particularly good coverage of our 
expected signal. 


There is a modest trend for points with smaller Yukawa 
couplings to decay later and therefore to have lower fre- 
quency peaks. For the signal to be observable, the Q- 
balls must decay when the temperature falls in the range 
20 GeV < Taec < 2 107 GeV. While the upper bound is 
frequently satisfied even for large reheating temperatures, 
a low reheating temperature is often preferred to avoid 
overproduction of gravitinos, although the exact bound 
on the reheating temperature depends on the mass of the 
gravitino M954]. We have imposed Tr < 10° GeV for 
all benchmarks. 


Conclusions: Affleck-Dine baryogenesis is consistent 
with a range of inflationary and dark matter models 
(see, e.g., ), as a high-scale phenomena, it is dif- 
ficult to directly test. We have shown that a broad 
class of Affleck-Dine models produce a detectable grav- 
itational wave signal within the range of the Einstein 
Telescope and/or DECIGO. Such signals are a conse- 
quence of the sudden end of an early matter-domination 
epoch, which occurs if the Q-balls from the fragmenta- 
tion of the Affleck-Dine condensate are sufficiently long- 
lived. A low reheating temperature, motivated by the 
gravitino problem, ensures a signal within the observable 
frequency range, but we find that this is not a require- 
ment. Thus, if a signal is observed, the cause is one of two 
known scenarios- an early period of Q-ball domination, 
which is a natural outcome of Affleck-Dine baryogenesis, 
or an early period of light primordial black hole domina- 
tion [14]. 


w (GeV) v (GeV) Yp 


r 


To (GeV) Ne Veft 


6.66 x 10° 3.80 x 107° 1.11 x 1078 
8.45 x 10° 1.92 x 10° 1.36 x 1078 
9.95 x 10° 7.21 x 10° 2.10 x 1078 


1.56 x 1078 4.59 x 10° 1000 0.024 
2.76 x 1077 8.04 x 10° 1000 1.4 x 107° 
1.38 x 107° 3.56 x 10° 1000 2.9 x 107° 


TABLE I. Parameters used in our three benchmark points in Fig.|1| In addition to the Q-ball parameters w and v, Ypg is the 
initial charge asymmetry after fragmentation which occurs at temperature To, No is the average initial number of Q-balls per 
Hubble volume after fragmentation, and r ~ Yp is the ratio of the asymmetric component. Note that the Yukawa couplings are 
equal to that of the Standard Model bottom quark, up quark, and electron in the top, middle, and bottom rows. Additionally, 


we have taken g» = 106 in our analysis. 


Qo Teq (GeV) Tace (GeV) 
1.14 x 107? 6.34 x 10° 1368 
1.47 x 10” 55520 138 
5.18 x 10° 20050 458 


TABLE II. Calculated quantities for the three benchmark 
points infi] Teq is the temperature of Q-ball-radiation equality 
and Taec is the temperature of Q-ball decay. 
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1. Explicit Potentials 


To keep the discussion as general as possible, we 
have avoided specifying a potential. In this appendix, 
we discuss Q-balls in both gauge-mediated and gravity- 
mediated scenarios. 

In the gauge-mediated supersymmetric scenario the 
potential is 


4 |e/? l aie 
V(®) = m* log (14+ + =|)’, (19) 


m2 


plus a small CP-violating term. When the second term is 


negligible, the resulting Q-balls are thick wall Q-ball [27], 
for which the analysis presented here is inapplicable. 
However, scaling arguments favor a Q-ball domination 
epoch. Fragmentation tends to produce one large Q-ball 
in each Hubble volume [20], resulting in a sharply peaked 
mass distribution at large masses, which tend to be long- 
lived. Furthermore, the radius now scales as Q!/4, and 
after accounting to the scalings of the VEV and energy 
per unit charge with Q, we expect the Q-ball decay rate 
to scale as Q'/4. This is suppressed compared to the 
thin-wall rate, and therefore, the Q-balls will tend to be 
longer-lived. T'Q-pan/H@ then scales as Q-3/ 4 which 
increases as the charge decreases, leading to the rapid 
matter-to-radiation transition. We plan to address this 
scenario more fully in future work. 

When the second term in is not negligible, then 
the resulting Q-balls are in the thin wall regime, although 
since the energy per charge w is independent of the charge 
Q, the equilibrium state is not one Q-ball per Hubble 
volume. 

Alternatively, one can consider gravity-mediated 
SUSY breaking, in which case the Affleck-Dine conden- 
sate has the potential 


2 

V(®) = m?|0|? (1 + Klog (=5)) + 1 j5j6, (20) 
where A is an effective scale for the higher-dimensional 
operator, m is the mass of the scalar field, and K ~ —0.01 
to —0.1 is a one-loop correction. Since K < 0, Q-balls 
can be formed, and in the thin wall limit, the VEV inside 
the Q-ball is [I] 


K 1/2 
vx (aum K) ; 


from which an expression for w can be found. 

In Table [IIT] we show one choice of parameters for the 
gauge-mediated potential (left) and gravity-mediated po- 
tential (right) for each benchmark set of parameters dis- 
cussed in the text. That is, in each row the Q-balls pro- 
duced have the same VEV v and energy-per-unit-charge 
w as in corresponding row in Table |I| For the gauge- 
mediated scenario, we have ensured that the sixth order 
term is relevant so that we are in the thin-wall regime. 
We see that in all cases the scale of the effective operator, 
A, is well above the reheating temperature. 


(21) 


2. Gravitational wave spectrum 


We here outline the gravitational wave spectrum from 
an instantaneous transition from matter to radiation 
domination, following Ref [3I]. We define the following 
functions: 


22 max 
80 (%, Umax) i (e | 2 


Bas — | + ==- v3) 


ax 


Gauge-Mediated Gravity-Mediated 
m (GeV) A (GeV) m (GeV) A (GeV) 
6.07 x 107 4.54 x 10” 2.16 x 10° 4.23 x 10'° 
1.76 x 10” 1.10 x 10+ 1.13 x 10° 2.07 x 1078 
3.61 x 107 1.35 x 1014 1.63 x 10° 2.01 x 1014 


TABLE III. Parameters for the potentials and which 
produce Q-balls corresponding to our three benchmark points. 
The first row corresponds to our first benchmark in Table 
the second row corresponds the second benchmark, and 
the third corresponds to the third. For the gravity-mediated 
scenario, we have fixed k = —0.05. 


Tmax Tmax 
Olzr-2 © |2 
| al | v3 
along with 


sie) = f poe oc) =- f n (8) 


zZ z 


The gravitational wave spectrum involves a resonant con- 
tribution, Q,.;, an infrared contribution, QyR, and a non- 
resonant UV contribution, Quy, such that the total spec- 
trum is given by 


Naw (z, Timax) — Qres+ 


3A2z8 40; [3]? + (T — 25;[$])? 22nax 1 2ns 
s™ max 217+2ns 625(3 + 2ns) r 
m 2(ns—1) 
x (=) (Oi 8 [max = x] Is Quy [x D Lmax]) , 


(24) 
where we define X = £/Zmax to write 


1 
(2 + ns) (3 + ns) (4+ ns) (5 + 2n,)(7 + 2ns) 


Or = 


x (1536 — 6144X + (7168 — 1920n, — 256n2) X? 
+(5760n, + 768n2) X? 

+(1328n, + 3056n? + 832n3 + 64n4)X* 

— (7168 + 12256n, + 7392n? + 1664n3 + 128n4)X* 
+ (7392 + 10992n, + 5784n? + 1248n3 + 96n2).X° 
— (2784 + 3904n, + 1960n? + 416n3 + 32n4) X7 

+ (370 + 503n, + 247n? + 52n3 + 4n4) XS 

— 256 (1 — X)® [(6 + 6(2 +n,)X 


+(2+ns)(5 + 2ns)X7] (1 - =x) | ’ (25) 


and 


Quv = 2(2— X)*T[4 + 2n,] 


X’ 4X? (2- X}? 24(2- X) 
j (z Fam) +m] ) p2) 
Q 


T[9 + 2n] 
2.34/33" TE 
res = ggg x 23ra? XY AS 0(2, Smax) 


1 2 
x C Tooo o ] 


3 S0 (x, Tmar) ] 


1 
32F Pea 
Fily nag 3 


3 5 Solz, max)” 
—so(T, max) 2 Fi [5 Ns, 2’ of 3 ) )) 
(27) 


where F} is the hypergeometric function. To convert to 
a frequency spectrum, simply take x, = 1/k,7, and 


Qaw (f)h? = 0.39h70,. 


T Te 
ees E x 10724 (ana) f, aa . (28) 


